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Abstract. In this paper we introduce a new PDE model in fre- 
quency space for the inertial energy cascade that reproduces the 
classical scaling laws of Kolmogorov's theory of turbulence. Our 
point of view is based upon studying the energy flux through a 
continuous range of scales rather than the discrete set of dyadic 
scales. The resulting model is a variant of Burgers equation on the 
half line with a boundary condition which represents a constant 
energy input at integral scales. The viscous dissipation is modeled 
via a damping term. We show existence of a unique stationary 
solution, both in the viscous and inviscid cases, which replicates 
the classical dissipation anomaly in the limit of vanishing viscosity. 

A survey of recent developments in the deterministic approach 
to the laws of turbulence, and in particular, to Onsager's conjecture 
is given. 



1. Motivation for the Model 

1.1. Onsager and Kolmogorov. The Euler equations for the motion 
of an incompressible, inviscid fluid are 

(1.1) — + (u ■ V)u = -Vp 

(1.2) V-u = 0, 

where u(x, t) is a divergence free velocity vector and p(x, t) is the inter- 
nal pressure. We consider the system in 3 spatial dimensions and we 
assume that the domain is either periodic or the entire space M 3 . To 
obtain the energy equation we multiply ( II. ip by u and integrate, using 
( P to give 

(1.3) il = - /(u-Vju.^. 
2 at I I 
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We define the total energy flux il by 

(1.4) il = (u- V)u ■ udx. 



For smooth solutions we can integrate by parts and use (11.21) to conclude 
that il = and hence energy conservation holds, i.e. 

(1.5) J \u{x,t)\ 2 dx = J \u(x,0)\ 2 dx for t > 0. 

However, in the context of turbulent flows in the limit of vanishing 
viscosity, it is appropriate to consider the Euler equations in the sense 
of distributions and impose only minimal assumptions on the regularity 
of the velocity field u. In the absence of sufficient smoothness we cannot 
integrate by parts in (11.41) or even make sense of (11.4j) and ensure that 
IT = 0. Conservation of energy might then be violated. Hence it is 
of interest to ask what are the minimal regularity assumptions on the 
velocity that ensures that (11.51) holds. 

Observing that the integrand in (11.41) is cubic in u and contains 
one spatial derivative suggests that if u has Holder continuity h > 
1/3, integration by parts is justified and n = 0. In fact this was 
the conjecture made many years ago by Onsager in his seminal paper 
on statistical fluid dynamics [17] . More precisely, he conjectured that 
(a) every weak solution to the Euler equation with smoothness h > 
1/3 conserves energy and (b) there exists a weak solution with h < 
1/3 which does not conserve energy. Such putative energy dissipation 
due to the irregularity of the flow is called anomalous or turbulent 
dissipation. A detailed historical account of Onsager's theory is given 
by Eyink and Sreenivasan [10] . 

All physical fluids are viscous, if only very weakly so. Turbulent 
fluids are believed to be described by the Navier-Stokes equations 

du 

(1.6) - + (m-V)« = -VP + i/A« + / 

(1.7) V-m = 0, 

where u, which could be very small, is the coefficient of viscosity and / 
is an external force which supplies energy into the system. The "classi- 
cal" Kolmogorov theory of turbulence predicts that energy dissipative 
solutions to the Euler equation may arise in the limit of vanishing vis- 
cosity for "generic" viscous flows that are governed by (II. 6p - (11. 7p . In 
homogeneous, isotropic turbulence the mean kinetic energy per unit 
mass is defined by E = ^(\u\ 2 ) while the energy density spectrum 
is defined by E(k) = \-^{\ u <k\ 2 )- Here u <K denotes the filtered ve- 
locity field containing all frequencies below a wave number k. Hence 
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E = J °° E{n)dK. The mean energy dissipation rate per unit mass is 
defined by 

(1.8) e v = (HVwf), 

where u u is a solution to ( II. 6p - (II. 7p . Kolmogorov [15] predicted 
that the energy cascade mechanism in fully developed 3-dimensional 
turbulence produces a striking phenomenon, namely the persistence of 
non-vanishing energy dissipation in the limit of vanishing viscosity, i.e. 

(1.9) \ime v ^e>U 

where e is the anomalous dissipation rate for the inviscid Euler system. 
The positivity of the limit in (II. 9p is referred to as the dissipation 
anomaly. Let us now assume that / = f <Kf (i.e. / has finite Fourier 
support) and that solutions to ( II .6p tend to a statistically stationary 
state with uniformly bounded mean energy. We multiply (jl.6p by u u <K 
and obtain 

(1.10) H^) = -^<|V<J 2 ) + </•<„). 

If k > Kj we have (/ • u v <K ) = (f ■ u v ) = e u > 0. On the other hand, 
by Bernstein's inequality, u(\ Vm"J 2 ) < uk 2 {\u v \ 2 ). Since the energy is 
uniformly bounded by assumption, we obtain from ( ll.lOp that 



(1.11) YmiUJn) = lime" = e. 

Thus in the limit of vanishing viscosity the average solution of the 
forced Euler equation inherits the anomalous dissipation rate e. 

As Frisch [12] describes, a self-similarity hypothesis on the velocity 
increments in small (spatial) scales implies that the energy spectrum 
as a function of wave number k has the power law 

(1.12) E(k) ~ £ 2//3 /t~ 5 / 3 

in the "inertial" range k G [/c/,ACd]. Here is the Kolmogorov dissi- 
pation wave number given by = (e/V 3 ) 1 / 4 and Kf = max{|fi;| : k G 
supp/}. This power law is known as the K41 turbulence model. Al- 
though the 5/3 power law is consistent with much physical data, there 
are also experiments which indicate turbulent regimes with alternative 
power laws. In fact, Kolmogorov's 1941 theory requires that the local 
velocity fluctuations are uniformly distributed over space. However, in 
reality dynamical stretching of the vortex filaments in 3-dimensional 
flows leaves some regions of the fluid domain with moderate turbulent 
activity and other regions with intense activity. This so called spatial 
intermittency should reasonably be accounted for in the description of 
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the scaling laws. The expressions for E(k) and that incorporate the 
dimension D of the effective dissipation region are 

(1.13) E{k) ~ e 2 /3 K -(8-D)/3 
and 

(1.14) K d ~ (e/^) 1 ^ 

where Z) £ [0, 3]. Thus the classical K41 model corresponds to D = 3, 
i.e. uniform distribution over 3 dimensional space, while D = cor- 
responds to a fully intermittent model where energy cascades through 
scales and dissipates only on points. 



1.2. Onsager's Conjecture and Besov Spaces. In the past few 
years there have been a number of articles that address part (a) of On- 
sager's conjecture. These include articles by Constantin et al jl], Eyink 
[H], Duchon and Robert [7J. It was shown that appropriate function 
spaces to examine the Euler equations in the context of Onsager's con- 
jecture are Besov spaces. In such spaces the notion of energy balance 
when the velocity is "a little smoother" than Holder h > 1/3 can be 
made precise. These are the natural spaces to work with in terms of 
a description of the energy flux phrased by a Littlewood-Paley decom- 
position which provides detailed information concerning the cascade of 
energy. Recently Cheskidov et al [2] obtained the largest Besov space 
where conservation of energy is ensured for the Euler equation. We 
note that to date there are no examples of Euler flows that possess 
some smoothness and confirm the second part of Onsager's conjecture, 
although there are examples of "very weak" Euler solutions that violate 
the energy balance condition [51 [HI [19] . 

We recall the definition of a weak solution of the Euler equation. 
A vector field u £ C^([0,T] : L 2 (IR 3 )) is a weak solution of the Euler 
equations with initial data u £ L 2 (IR 3 ) if for every compactly sup- 
ported test function ip £ C£°([0, T] x IR 3 ) with V x -ip = and for every 
< t < T, we have 

(1.15) u(t)ij){t) - u{0)i/>(0) - I u-d s ipds= [ [ u-V^-udxds 

Jo Jo Jr s 

and V x . • u = in the sense of distributions. 

We define the Littlewood-Paley energy flux EL, through a sphere in 
frequency space of radius 2 J as follows. For any divergence free vector 
field u £ L 2 (R 3 ) we define 



(1.16) 



S J u = u*F- 1 {4){-2- ] )) 
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where ^(O is a smooth nonnegative function supported in the ball of 
radius one centered at the origin and such that = 1 for £ < 1/2 
and J 7 is the Fourier transform. We then define Uj as 



(1.17) Uj = - u-S/SjU-udx. 

Using the test function S 2 u in the weak formulation of the Euler equa- 
tions we obtain 

(1.18) ~4-l|s>||| = -il 

V J 2dt \\ J 112 J 

Cheskidov et al pj prove that the Littlewood-Paley energy flux of a 
divergence free vector field u G L 2 satisfies the following estimate: 



//•ll 3 

' i II 3 



(1.19) |n,| < 

i=-l 

where Uj is the j-th Littlewood-Paley piece of u defined by 

Uj = Sj + \u — SjU. 

It follows from (1.19) that 

(1.20) limsup|rLj < Climsup(2 j ||u i ||^). 

Furthermore, an important feature of the bound (1.19) is that it is 
quasi- local in the sense of rapid decay when \j — i\ is large. 

1/3 

We define the Besov space B 3 ' CQ to be the space of all tempered 
distributions «6l 3 for which 

(1.21) limsup2 i/3 ||u j || 3 = 0. 

Hence, from (1.18) and (1.20) we obtain the following result. 

Every weak solution u to the Euler equation on a time interval [0, T] 
which satisfies 

(1.22) lim / 2 i ||w j (t)|| 3 ; = 

conserves energy on the entire interval [0, T\. In particular, energy is 

1 /3 

conserved for every solution in the class L 3 [0, T]; B 3co )f]C w ([0, T\\ L 2 ). 

In order to see more transparently the connection between (1.22) and 
the smoothness 1/3 predicted by Onsager we rewrite (1.22) as follows: 

(1.23) lim -r^r / f \u(x) -u(x - y)\ 3 dxdt = 0. 

\y\ Jo Jr3 
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Hence the solution to the Euler equation needs to be a little better than 
1/3 Holder continuous in the space-time average to ensure that energy 

1/3 

is conserved. We call the Besov space B 3 oo Onsager critical. This is the 
space which contains distributions u G M 3 where lim^oo 2 j, / 3 ||m j || 3 is 
finite, but not necessarily zero. This is a critical space in which energy 
conservation for the Euler equation might be violated. 

Applying the bound on the energy flux given by (1.19) to the Navier- 
Stokes equations gives a sufficient condition for the energy equality to 
hold, namely. 

Let u v G C w ([0,T];L 2 (R 3 )) n f 2 ([0, T]; H^R 3 )) be a weak solution 
to (1.6)- (1.7) with 

(1.24) lim / 
Then u v satisfies the energy equality 

(1.25) \\u u (t)\\l + 2u J \\Vu u (s)\\ 2 2 ds = \\u u (0)\\ 2 2 + 2 \ f-u u (s)ds. 
In particular, (1.25) holds if u G L 3 ([0,T]; if 5 / 6 ). 



2 3 '||u^(t)||ldt = 0. 

'o 







1.3. Littlewood-Paley Framework for Intermittency. Let u v be 

a Leray-Hopf weak solution to the Navier-Stokes equations (11.61) - (II. 1\ . 
We denote by (•) the long time average. We define the Littlewood-Paley 
energy spectrum of u v by 

(IKII2) 

(1.26) E lp {k) = 

K 

for frequencies k G [2 j , 2 j+1 ] and we define the mean energy dissipation 
rate by 

(1.27) e" = i/<|| V«12). 

If a family of individual realizations {u u } u<uo verifies Kolmogorov's 
hypothesis that e v — > e > 0, then the locality of the flux which is 
exhibited in the bound (1.19) suggests the following 

(1.28) 2*<||ti;||jj)~ e 

for all j sufficiently large. Here u° = lim^^o^^- 111 other words, the 
limiting solution u° to the Euler equation is "on average" in the Onsager 
critical space -63 ^. 

Eyink [9] showed that B\^ is consistent with the multi-fractal in- 
termittency models of Frisch and Parisi [11]. Within the Littlewood- 
Paley framework we can model the intermittency correction, (see (11.131) 
and (11.141) ) by assuming the relationship between e and ||w?||3 given in 
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(I1.28P and fixing the saturation level in Bernstein's inequalities. More 
precisely, in 3 dimensions we have 

(1-29) lk-||3<2 i/2 ||%|| 2 . 

Assuming that the region of active turbulence is bounded, (say, on a 
torus) we also have 

(i- 30 ) Ikih ^ \\uj\\ 3 - 

Hence for 2 J £ [tzf, Kd] we write 

(1.31) K-lla ~ 2^11^112 



for some c in the interval [0, 1/2]. Then from (ll.26p . ( 1 1 . 2 8 [) and fll.3ip 

we recover the energy spectrum law 

2/3 

(1-32) E lp (k) ~ JL^ 

with 2 J being identified with k. 

The analogy between f 1 1 . 1 3 j) and fll.3ip requires that D = 3 — 6c. So, 
the fully saturated Bernstein's inequality (i.e., c = 1/2) corresponds 
to a uniform distribution of modes Uj in each dyadic shell and hence 
strong localization in space (i.e. D = 0). On the other extreme, c = 
corresponds to a uniform distribution of u in physical space space (i.e. 
D = 3) and localization in frequency space which corresponds to the 
classical K41 model. 



2. A Continuous Model for the Energy Flux 

Although there is abundant empirical evidence for Kolmogorov's hy- 
pothesis that \im u ^ Q e u e > 0, this has not been rigorously proved 
for the Navier-Stokes to Euler limit. It is therefore of interest to ex- 
amine simpler models that retain some of the essential features of the 
fluid equations and yet are tractable enough to allow a proof of Kol- 
mogorov's hypothesis. We now propose a PDE model for the turbulent 
energy spectrum in frequency space. We choose the scaling to include 
the intermittency correction that we described in sections 1.2 and 1.3. 
To motivate the model we start with a fully local version of the flux 
given by the bound in the inequality (1.20), namely 

(2.1) n,- ~ a'lMl. 

We further assume that c in f ll.3ip is independent of j. We thus obtain 

(2.2) Ilj ~2 {3c+1)j \\u J \\l 
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We now make a further step from the discrete expression for flux ( 12.2ft 
to a continuous one by looking at the energy density function a(n, t) 
defined by 

\\ U (t)\\2 = / 2 (^, t)dK. 







Assuming that a "does not vary much" in each dyadic shell, or disre- 
garding energy density variations in each dyadic shell, we obtain 



ll M jll2 ~ / a 2 (n,t)dn ~ Ka 2 (n,t). 

for all k G [2 j ,2 j+1 \. Thus, (E2D becomes 

(2.3) n(/c) = K 3c+ fa 3 (/«,t). 

Going back to (11. 101) we assume that = 1, and the energy in the 
sub-inertial range in negligible. So, with ( 12. 3 j) at hand we can write 
the energy balance relation as follows 

^Qf a 2 (£,t)d/) = -/t 3c+ fa 3 (/t,t) -z/jT £ 2 a 2 (e,t)di. 

Differentiating in k and cancelling a on both sides we obtain the fol- 
lowing PDE 

5 3 5 

(2.4) a t = -(- + 3c)/t2 +3c a 2 - 3^ +3c aa K - vn 2 a 

We supplement this equation with the boundary condition 

(2.5) a(l,t) = e 1/3 . 

Here e represents the energy input rate coming from an external force. 
We thus disregard any particular detail of energy production and sim- 
ply model it with our boundary condition. The input rate e will sub- 
sequently be shown equal to the energy dissipation rate, hence the no- 
tation. Equation ( 12 .4p can be easily simplified by rescaling the energy 

5 

density a to b = n& +c a. Thus, ( 12. 4p becomes the following equation 

(2.6) b t = -3n a bb k - isK 2 b, k > 0, 

where a = | + 2c. The appropriate range of a is [|, |], which exactly 
corresponds to the classical range of the energy density power laws with 
the spatial intermittency correction. 

Remark 2.1. Equation ( 12. 6 p is a variant of the much studied Burgers 
equation. The one dimensional Burgers equation can be viewed as the 
most basic nonlinear PDE that has the bilinear structure of the non- 
linearity of the Euler and Navier-Stokes equations. It can be invoked 
as a model for one dimensional compressible fluids. However there is 
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no clear physical basis for using Burgers equation in physical space as 
a model for turbulence. On the other hand, as we have argued, the 
locality of the energy flux manifested using Littlewood-Paley theory, 
motivates ( 12. 6 j) as a PDE model for the turbulent cascade in frequency 
space. 

Remark 2.2. In the past few decades a number of "toy models" for 
turbulence have been studied to test Kolmogorov's theory. In par- 
ticular, the derivation of the classical Desnyanskiy-Novikov discrete 
model, [B] , follows a similar path. The flux there is modeled by taking 
Uj = 2 d iaj(ij + i, where a| represents the total energy in the j-th dyadic 
shell, while d is an inter mittency parameter with the appropriate range 
of values. The model is thus an infinite system of ODEs given by 

(2.7) ja 3 + u2 7j a j - 2^-^aJ_ 1 + 2 dj aj a j+1 = 0, j = 0, 1, 2, . . . , 

where a_i = 0. This model, as well as its inviscid versions, has been 
extensively studied by Katz and Pavlovic [13], Cheskidov and Fried- 
lander [3], Kiselev and Zlatos [H], and others. To our knowledge the 
PDE model we present is the first continuous model. 

In Section [3] we examine the inviscid [y = 0) form of (12. 4p with 
boundary condition (12. 5p . We prove that there is a unique fixed point 
which is a global attractor. Moreover every solution reaches it in finite 
time. The inviscid equation exhibits anomalous dissipation and the 
average energy spectrum has the power law e 2 ^ 3 K~ a . 

In Section H] we turn to the viscous model v > (12 .4p . which is 
in essence a damped Burgers equation. Again the PDE has a unique 
fixed point and this converges to the inviscid fixed point as v — > 0. The 
viscous fixed point reproduces Kolmogorov's energy density spectrum 
in the inertial range and it becomes zero identically after the dissipa- 
tion wave number k^. We further consider the Leray regularization of 
equation (12 .4p . We show that all bounded solutions of the regularized 
equation converge pointwise and in the metric of L 2 -space to a fixed 
point which in turn converges to the fixed point of equation (12 .4p in the 
limit as the regularization parameter goes to zero. The average dissipa- 
tion rate for the viscous system converges to the anomalous dissipation 
rate for the inviscid system giving an example of the dissipation anom- 
aly predicted by Kolmogorov. 
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3. Inviscid case 
In this section we study the inviscid version of the model ( 12 Ah 



(3.1; 



a 

a t = — 3aK a -— (ft2 a), k > 1, 

OK 



a(l,t) = e 1/3 , 
a(n, 0) = a (n) > 0, 

where a G [5/3,8/3]. Note that the energy equality 

d 1 f°° 

(3.2) -- J a(K,t) 2 dK = e 

is satisfied on some interval t e (0, T) provided the solution satisfies 
the following smoothness condition (cf. dl .22[) ) : 

r T 

(3.3) lim / n 3a/2 a 3 dt = 0, 

which is an analog to ( ll.22p . So, e > represents the energy input rate 
in this model. The unique fixed point of ( 13. ip is given by 

(3.4) A\k) = e 1 / 3 ^ 2 

We note that the fixed point does not satisfy ( I3.3p . Moreover, it does 
not satisfy the energy equality since 

(3.5) ±-J A°(K) 2 dK = 0^e. 

The anomalous energy dissipation rate is the difference between the 
energy input rate and the time derivative of the total energy. Thus 
using (3.2) and (3.5) we observe that the anomalous dissipation rate 
for the fixed point is exactly the energy input rate e. We will show 
that this also holds for every other solution asymptotically in time. In 
order to do this we will prove that A is a global attractor. 
We use the following change of variables: 

£ = t) = e~ 1/3 C^a(r\ K 1/3 7*), 

where 7 = -^r- Then ( 13. ID reduces to Burgers equation 

(3.6) w t = wwe, < f < 1 

with w(l,t) = 1, 0) = u>o(£) > 0. We extend it to Burgers 
equation on the whole real line 

(w t = wwz, i E R, 
1 } \w(Z,0) = wo(0, 
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where 

mo = I mo, o < <e < i, 

The weak solution to (13. 7p is expressed using the Lax-Oleinik for- 
mula. Let 

h(y)= [ V w (0^. 
Jo 

For all t and for all but at most countably many ^Gl, there exists a 
unique ?/*(£,£), such that 

min U(y)} = f(y*)- 

y 

where 

f(y) = ^f^-h(y). 

Then 

y.{Z,t)-S 



(3.8) w(£,t) 



t 



is the weak solution to flU]). Given £ G [0, 1], let t > 2(1 — £). We 
will show that y*{£,t) = £ + t and consequently w(£,t) = 1. Indeed, 
let y = £ + s. First, consider the interval s > 1 — £. Then 



s 2 ,1 



f(y) = Yt~ J Q M0dt-z-s + i> f(t + t) 

provided s ^ t. Since t > 2(1 — £), it follows that 

/(£ + *)<- fw {i)di, 

Jo 

and hence on the interval s < 1 — £ we also have 

f(y) >Y t - f MO d ^~J MO d£ > m + *)■ 

Therefore, (13. 8p implies that w((,, t) — 1 for £ G (0, 1), t > 2 and hence, 
returning to the original variables, 

(3.9) a(«,t)=A°(K), t>|^ 1/3 7. 

Hence the average energy spectrum for solutions to the inviscid model 
is about e 2 / 3 K~ a . 
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4. Viscous case 
In this section we study the viscous model (12 Ah 



(4.i; 



d , at v q 

a t = —3an —{K 2 a) — vk a, k > 1, 

OK 

a(l,t) = e 1/3 , 

a(n, 0) = ci (k), a o ^ -^ 2 ; a o > 0, 



where a G [5/3,8/3]. There exists a unique fixed point to (14.1 ft given 
by 



(4.2) A"(k) 



-a/2 



6^ + 



V 



K 



3-a\ 



1 < K < K d , 



3(3 - a) v 

0, K> Kd, 

where k^ is Kolmogorov's dissipation wavenumber described in Section 
1.1. For the model it is explicit and given by 



(4.3) 



1 + 



3(3 — a)e 



1/3 



i 

3-q 



To see the parallel with the classical expressions for Kd we note that 
for v small one has 

5 
3' 



Kd 



e \ 4 



v 



for 



a 



and 



r, for 



a = -. 
3 



In the limit of vanishing viscosity we immediately obtain from (4.2) 
and (3.4) that 

(4.4) A v {k) -> A°(k), 

uniformly on any finite interval [1,acq] as v — > 0. With a little more 
effort we can show that the convergence also takes place in L 2 ([l, oo)). 
Indeed, 



\A U {K) - A {K)\ 2 dK < 6 2/3 / K~ a dK+ 

J K d 



K 



d,K. 



9(3 -a) 2 J 1 

Since Kd — > oo we see that the first integral vanishes as v — > 0. The 
second integral behaves like v 2 for a < 7/3, like z/ 2 log(^) for a = 7/3, 
and like z/ 0-1 ^ 3-0 ) for 7/3 < a < 8/3. So, within our range of a the 
second integral vanishes too as v — > 0. 
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In order to study the time dependent solutions to the viscous system 
(14.1 p we utilize the same change of variables as in the previous section 

t = t) = e-WC^aiC, K V3 7*), 

where 7 = -4r- Then (14. ip reduces to the damped Burgers equation 

(4.5) w t = wwf: - ^C 2l w, < £ < 1 

with w(l,t) = 1, w(£,0) = w (0 > 0, fiew*(£)d£ < 00. Here 



jj, = \ve 1//3 7 . The equation has a unique fixed point 



(4.6) w(o = { 2 7 -i l 

0, < z < 



where 
(4.7) 



2 7 - 1 
1 + — 



1-27 



Note that ^ — >• as n — >• 0. Note that «d = £J 7 . 

Since there are discontinuous solutions to the damped Burgers equa- 
tion ( I4.5p . we now consider a Leray-type regularization of the equation. 
Such regularizations have been used to approximate weak solutions to 
the Burgers equations [H HE] . Consider the following regularized equa- 
tion: 

(4.8) w t = vsw^ — fi£~ 2j w, vg = w * 4>Si 0<£<1,£>0, 

with the boundary conditions w(l,t) = 1, w(0,t) = 0. Here <f>g(£) = 
-g4>(£/8), where </>(£) is smooth, nonnegative, and such that f </>(£) d£ = 
1, supp0 = (—1,0). First, note that there exists a unique fixed point 
w = Ws(£) to ( 14. 8p . which is smooth, nonnegative, monotonically in- 
creasing, and with H^O) = 0. Now consider characteristics r)s(t), which 
are solutions to 

(4.9) - V6 (t) = -v 6 (r 1s (t),t). 

It is easy to see that characteristics do not intersect on (0, 1). Indeed, 
take e G (0, 1) and consider two solutions r)' s (t), r)g(t) to ( 14.9D . As long 
as they belong to (e, 1) we have 

j t \m -m\ < s- 3/2 \m -m\ (^(e, W 1 

1 /2 
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Since £ 2 w 2 (t;,t) d£ is non-increasing, the characteristics do not in- 
tersect. 

Along characteristics we have 

^w s (rjg(t),t) = -fMrjg 2j ws(rjs(t),t). 

Moreover, along every characteristic curve that starts from the bound- 
ary £ = 1 we have that Wg is equal to the fixed point Wg, i.e., 

w s (m(t),t) = W s (r]s(t)), 

provided rjg(to) = 1 for some to > 0. Now consider the characteristic 
curve rfg{t) with rjg(0) = 1. Note that rfg (t) is decreasing, positive, 
r]g(t) — > as t — )■ 00, and 

(4.10) w s {Z,t) = W s {0, 

provided £ > rjg(^,t). Note also that rjg (t) < T)§ 2 (t) provided 5\ < 62- 
Then consider the particular characteristic curve 77° (t) corresponding 
to 5 = 1. Let T(f) be such that rji( T (0) = £• Let £0 G (0, 1). Then 
for all 5 < 1, f <E [f , 1], t > T(£ ) we have 

(4.11) = W*(0- 

So the fixed point of the regularize equation Wg is a global attractor. 
Furthermore, as 5 — > the fixed point Wg converges to the fixed point 
W given by (14. 6p . 

We now return to equation (14. ip for a(£, t). We will study the energy 
dissipation in the limit of vanishing viscosity. For convenience, solu- 
tions to the model (14. ip with viscosity v > will be denoted by a"(t). 
The fixed point (which is unique in both viscous and inviscid cases) 
will be denoted by A u . We assume that a u (-,0) e L 2 and a v (t) ->■ A u 
in L 2 . We denote the energy and enstrophy norms of a u by 



00 

\a v {t)\:= J a v {n,tfdK, := / K 2 a v {K,tf da. 

Due to the energy inequality we have 

^K(T)| 2 -^K(0)| 2 <-zJ, / \\a"(t)\\ 2 dt + e. 



lim sup - / u\\a v (t)fdt < e. 



Hence, 



On the other hand, note that the fixed point A v satisfies the energy 
equality 

v\\A v t = e. 



A CONTINUOUS MODEL FOR TURBULENT ENERGY CASCADE 
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Now for any r] > 0, there exists N, such that 

r»jV 



K 2 A u (n) 2 dK > u\\A u \\ 2 -rj. 



Since a u (t) — > A v in L 2 , we have 

fT r N rN 



i pi piv pj\ 

liminf-y J VK 2 a u (K,t) 2 dKdt>v J K 2 A u (K) 2 dK>iy\\A u \\ 2 -r]. 
Therefore, 

liminf- I v\\a v {t)\\ 2 dt > v\\A v \\ 2 = e. 

T^oo T J 



Then we obtain 



:=lim lim i / f dt > 0. 



i.e., in the limit of vanishing viscosity the energy dissipation rate is 
positive and equal to the energy input rate or anomalous dissipation 
rate for the inviscid model. 
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